Qubit Mediated Time Robust Entangling of Oscillators in Thermal Environments 
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We consider two separated oscillators initially in equilibrium and continuously interacting with 
thermal environments, and propose a way to entangle them using a mediating qubit. An appropriate 
interaction allows for an analytic treatment of the open system, removes the necessity of fine-tuning 
interaction times and results in a high tolerance of the entanglement to finite temperature. The 
entanglement thus produced between the oscillators can be verified either through a Bell inequality 
relying on oscillator parity measurements or through conditional extraction of the entanglement on 
to two mutually non-interacting qubits. The latter process also shows that the generated mixed 
entangled state of the oscillators is an useful resource for entangling qubits. By allowing for in- 
fluences from environments, taking feasible qubit-oscillator interactions and measurement settings, 
this scheme should be implementable in a variety of experimental setups. The method presented 
for the solution of the master equation can also be adapted to a variety of problems involving the 
same form of qubit-oscillator interaction 



INTRODUCTION 

Two of the simplest quantum systems are the qubit 
and the harmonic oscillator (oscillator for brevity). Their 
controlled interaction is realizable in a multitude of 
settings. These include cavity-quantum electrodynam- 
ics (QED), where atomic qubits interact with electro- 
magnetic fields circuit- QED, where superconducting 
qubits interact with microwave resonators Q or quan- 
tum nano-mechanics where mechanical oscillators inter- 
act with superconducting qubits [|[ and atomic qubits 
[H . It is generally perceived that the oscillators will have 
to be first cooled to a known pure state and isolated suf- 
ficiently from their environments before quantum effects 
in them can be seen. However, one can ask the following 
question in order to lower the demands for controlling 
quantum systems. How far can cooling and isolation be 
relaxed? We answer this question in the context of en- 
tangling two oscillators by a mediating qubit. This is 
important for those qubits which do not interact with 
each other. The entanglement of the oscillators could 
then be a resource for entangling qubits and their tele- 
portation [1,@,[7, 8]. Additionally, entangling mesoscopic 
oscillators, as available in quantum nanomechanics, fun- 
damentally probes the boundary of the quantum world. 

The most natural initial state of an oscillator is a 
thermal state. Moreover an oscillator continuously in- 
teracts with its thermal environment even while inter- 
acting with a qubit. Can one use the available qubit- 
oscillator interactions to entangle two oscillators in the 
above situation in a verifiable way? Surprisingly, though 
work on entangling oscillators initially in "pure states" 
is aplenty 0, HJt], H, G3 , only very recently the first ma- 
jor step in answering the above question has been taken. 
Refs.[ll|, ll2| have shown how to entangle two thermal 
state oscillators by a mediator in a pure state with some 
heuristic/approximate considerations of dccoherence in- 
duced by zero temperature baths (as opposed to a ther- 
mal environment). However, wherever thermal states are 
relevant, the thermal environment is also unavoidable. 



An exact analytic treatment of this problem by solving 
the appropriate master equation for a thermal environ- 
ment is required to go beyond heuristics and ascertaining 
how the openness of the oscillators and their bath tem- 
perature constrain an entanglement generation protocol 
and its outcome. In this paper, we show that the above 
analytic treatment can be achieved by choosing a specific 
form of the qubit-oscillator interaction which is naturally 
available in quantum nanomechanics and can be imple- 
mented with minimal effort in cavity/circuit-QED. To- 
gether with the inclusion of a thermal environment for 
the oscillator, our treatment also enables the analytic in- 
clusion of the dephasing of qubits that may be important 
in a quantum nanomechanical setting. The scheme we 
are considering for the entanglement preparation, rely- 
ing on the strength of resonant interactions, can be faster 
than that with cross-Kerr-like interactions [ll[ . We also 
find an allowed time window for the entanglement to ap- 
pear which becomes broad when decoherence and tem- 
perature are low. This can make the scheme more robust 
to errors in interaction times than most other schemes 
d, 0] . In addition, we show the possibility of reciproca- 
tion of the entanglement of the oscillators to two qubits. 
Thus the generated mixed state of the oscillator can serve 
as a resource for qubit based quantum information pro- 
cessing (QIP). 

Our results should be particularly useful for quantum 
nanomechanics, where oscillators are yet to be cooled to 
the ground state. In technologies wherepreparation of a 
ground state oscillator is available p], [J, using thermal 
states will surely lower the demands on cooling and iso- 
lation. Our scheme will show that oscillators in almost 
classical states can indeed be entangled to provide a QIP 
resource even when they are continuously (but weakly) 
coupled to a thermal environment. 

In addition to this, we believe the solution method 
presented here to be relevant on its own, as it could be 
applied to a variety of problems involving the same in- 
teraction Hamiltonian and the same decoherence model. 
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PROBLEM DESCRIPTION 

We start by presenting the schematics of our entangle- 
ment generation protocol without specifying the Hamil- 
tonian. At first, we let two identical oscillators interact 
with the same qubit, which we call the entangling qubit, 
either sequentially (such as for an atomic qubit flying 
through a pair of cavities) or simultaneously (such as for 
a superconducting qubit). In the former case we will 
assume the two interaction times to be the same. The 
states of the qubit are {\g), |e)}, while the oscillators are 
described by the annihilation operators a, b. Let the os- 
cillators be initially in a thermal state and the qubit in 
\g) so that the initial state of the complete system is: 

e -/3a t a e -,3b t b 

p(0) = \g){g\®p ai = \g)(g\® — 2 , (1) 

where Z = Tr{e-^ ata } = Tr{ e -' 3f,tb } and (3 = uj/{k B T). 
Throughout the paper, the unit of % = 1 is assumed. 
After an interaction between the qubit and the oscilla- 
tors for a time t the system will evolve to p(t). At this 
point, we measure the internal state of the qubit to ob- 
tain two possible outcomes f = g,e. As a result the 
oscillators are projected in the mixed state Pf(t). Our 
aim is to show that for an appropriate interaction, such 
a state may be entangled by an amount which depends 
on the temperature T and the strength of decoherence. 
One could try to detect the entanglement contained in 
the state pf, for example by verifying a Bell inequality 
violation by local measurements on the oscillators (we 
will discuss how to do this for our specific model). Al- 
ternatively, one could extract part of the entanglement 
on to an additional pair of qubits A and B (entangle- 
ment "reciprocation"), initially in the separable state 
Iff; ff) — \q)a ® \g)s by making mode a interact with 
qubit A and mode b with qubit B respectively. The re- 
sulting operation alone may not suffice to entangle A and 
B, as they will, in general, remain entangled with modes 
a and b. So we perform additional measurements of mo- 
menta p a — —i(a — a^); pb — —i(b — b^) or the pari- 
ties II a = (-l) a+a ; fl b = (-l) fet6 of the two oscillators. 
Note that no entanglement is created during such recip- 
rocation since the two parts of the system ( A, a and -B, b) 
are not directly interacting and the measurement is per- 
formed locally. Thus measuring A and B to be entangled 
is a sufficient condition for the state p/(t) to be entan- 
gled. It also demonstrates that the oscillators can store 
entanglement to be later extracted by qubits. 



SOLUTION METHODS 

We present the form of the qubit-oscillator interac- 
tion we use, first showing the case of a single oscillator 
and then generalizing to the case of two oscillators. The 
qubit and the oscillator are resonant at frequency uj. The 



Hamiltonian in the interaction picture is: 

H = a 1 (a + J), (2) 

where we have introduced the Pauli operators o~\ = 
\e){g\ + \g)(e\;a 2 = i(\e){g\ - |fl)(e|);<7 3 = |<7><<?| - |e)(e|. 
This Hamiltonian can arise in a Jaynes-Cummings sys- 
tem by driving it with an external laser (this revives 
the counter-rotating terms while preserving the strength 
of the coupling [loj ) . It also arises naturally in nanome- 
chanical oscillators coupled to charge qubits (although 
in the latter rotation of the qubit basis such that 

03 — * o~\ is needed). The time evolution operator is 

U(t) = e- lHt = D(-iait) (3) 

where D(a) = e aa '- Q a is the displacement operator 
[HI (we can use o\ inside this operator as [a\a, aia^] — 
[a, aJ] = 1). If the qubit is interacting with two oscillator 
modes a, b the previous expression is modified into: 

U(t) = Dai-ia^Dbi-iatf), (4) 

where Di is the displacement operator for mode i [l3j . If 
the entangling qubit is measured in \g) or |e), the result- 
ing operators acting on the oscillators (combining evolu- 
tion and measurement) are 

(g\U(t)\g) = \{D a {it)D b {it) + D a {-it)D b {-it)), 
(e\U{t)\g) = -\(D a (it)D b (it) - D a (-it)D b (-it)). (5) 

The above operators have increasing entangling capabil- 
ities as t gets larger, as can be seen for example by ap- 
plying one of them to a product of two coherent states. 
At this point we include the effect of a thermal environ- 
ment in our treatment, for which we first get back to 
the interaction of a single oscillator with a single qubit. 
Introducing decoherence, the dynamics can be described 
by the master equation: 

^ = -i[a 1 ( a + a^,p\+L(T)p + L i) p, (6) 

where L(T) is the Lindblad operator for the oscillator in 
a thermal bath at temperature T, which will, for cases 
we consider, be of the amplitude damping form (l3j : 

L(T)p = |(n(T) + l){2apa t - a ] ap - pa i a) + 

+ —n(T)(2a i pa — aa^ p — paa^), (7) 

where n(T) = (exp(j3) — 1) _1 and k is the damping rate. 
The operator L^p — ^(a\po~\ — p) will be present for the 
case of charge qubits whose dephasing can be significant 
(here o\ instead of 03 has been used because for a charge 
qubit the basis is rotated, as mentioned earlier). The 
model is analytically solvable, and we think it is worth to 
present here the solution method in full detail, as it may 
be applied to other problems involving the same form 
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of interaction Hamiltonian. We switch to the Wigner 
representation of the oscillator defined as 



W(a,t) = -jL J d- 



r\ e 



r/ a — r/a 



Tr{D(r,)p(t)}. (8) 



This representation is widely used in quantum optics, 
where the density matrix of a single (or multi) mode 
field is encoded into the real- valued Wigner function [l3[ . 
However, since in this paper we are dealing with a com- 
posite system (oscillator + qubit), the state of the system 
can not be represented by a single function; it is indeed 
encoded in a 2 x 2 matrix of functions, which we may 
call the Wigner matrix. The Wigner matrix component 
Wi j is defined as the Wigner function for the operator 
(i\p{t)\j) where i,j = e,g. Note that Wij is not neces- 
sarily real valued for i ^ j; nevertheless if p is a genuine 
density matrix then the condition Wij = W*^ holds, i.e. 
the Wigner matrix is hermitian for any (a,i). The mas- 
ter equation for p can be converted into a differential 
equation for W by using the following correspondences 
(for a rigorous derivation see for example [HI): 

ap^ (a+^d a ,jW; pa (a - ^d a ^j W; (9) 

a f P -> U* - ^dA W; pa) -> (a* + ^ W. 

(10) 

It follows that the Wigner matrix obeys the differential 
equation 

8 t W = -i{a + a*)[a u W] - l -{d a * - d a ){a 1 ,W}+ 



L{T)W + L^W, 



(11) 



where we introduced the anticommutator {01,02} = 
01O2 + 02O1. The Lindblad operator in terms of the com- 
plex variable a is 

L(T)W= ~ (d a a + d a *a* +2A(T)d a d a *)W, (12) 

where A(T) = n(T) + 1/2. Note that dUD is indeed 
a system of coupled differential equations for the four 
components of the matrix W, which can be separated by 
choosing an appropriate decomposition of W in the Pauli 
basis {1, <7i, (72, £73}. In fact we can see that the commu- 
tator [o"i, •] in (jlip is non- vanishing only if it is applied 
to the Pauli operators 02 , cr 3 , while the anticommutator 
{o"i, •} is non- vanishing only when it is applied to the op- 
erators 1,(71. It is then a simple step further to verify 
that 



[0-1,0-2 ± ia 3 ] = ±((7 2 ± ia 3 ) 

{(7i,l±(7i} = ±(l±(7i). 



(13) 
(14) 



We have at this point all the tools we need to write down 
a decomposition of the Wigner matrix in " normal modes " 



W = \ (wi(l + (7i )+ 1-2(1-0-1 )+ 
+ ^3(0-2 - io- 3 ) + ^4(0-2 + 10-3)), 



(15) 



which then yields four decoupled differential equations: 



d t vi 



:(d a * — d a )v\ + L(T)vi, 



(16) 
(17) 



d t v 2 = -^(da* - d a )v 2 + L(T)v 2 , 

d t v 3 = 2i(a + a*)v 3 + L(T)v 3 -jv 3 , (18) 
d t v 4 = -2i(a + a*)v 4 + L(T)v 4 - jv 4 - (19) 

We suppose that initially the oscillator is at thermal equi- 
librium with its environment at temperature T, while the 
qubit is in the ground state. In terms of initial conditions 
this translates into 

vi(0) - v 2 (0) = -*«a(0) - iv 4 (0) = W T {a), (20) 

1 Ic| 2 



Wr(a) = — — e A 
7rA 



(21) 



These initial conditions allow for a compact solution 
which we are going to present in full, while the solution 
for generic initial conditions can be expressed in terms of 
convolution integrals as discussed in Appendix. Equation 
(jTCJ) can be solved by choosing an ansatz of the form 



vi(a,t) = W T {a + X(t)). 



(22) 



We know that W T {a) is such that L(T)W T (a) = 0, and 
clearly d t W T {a) = 0. Thus, if @ holds, it follows that 

L(T)v 1 = ~(\d a + \*d a ,)v 1 , (23) 
d t Vi = (xd a +X'd e r)v v (24) 

By substituting these into Eq. (fl6|) and collecting terms 
proportional to d a v\, we obtain the differential equation 



2 2 ' 



(25) 



with initial condition A(0) = as a consequence of (|20p . 
Such equation can be easily solved, yielding 



A(t) = - (1 



(26) 



Following the same procedure we obtain 

v 2 (a,t) = W T (a- A(t)). (27) 
For ^3 we consider the following ansatz: 

v 3 (a, t) = iW T (a)e^ t)(a+a * )+iyW . (28) 
Similar considerations as before allow us to express 

L(T)v 3 = ^(-(a + a*)fi + 2A^ 2 )v 3 , (29) 



d t v 3 = ((a + a*)fi + v) v 3 . 



(30) 



This time, after substituting such expressions in Eq. (|18p . 
we can collect separately the terms proportional to av 3 
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and those proportional to V3, to obtain the pair of coupled 
differential equations 



(i = 2i- -fi, 
i> = kA/i 2 — 7. 



(31) 
(32) 



From (|2"0"j) we get the initial conditions /z(0) = u(0) = 0, 
then the solutions of Eqs. (|3"Tj) and (|32[) are 



fx(t) = 4A(t), 
v(t) = -yt + kA I h(t) 2 <1t 



(33) 
(34) 



Again, U4 can be solved similarly, yielding 

Vi(a, t) = -iW T (a)e~^ t){a+a ' )+u{t K (35) 

At this point we can go back to our original problem 
and generalize the procedure for a qubit interacting with 
two oscillators. We have in this case a master equation 
for the two-mode Wigner matrix W(a,f3). If we sup- 
pose that the two oscillators have the same temperature 
and the same value of re, the solution is obtained from 
the single mode case by replicating each function, i.e. 
Vi(a,t) — > Vi(a,t)vi(P,t). Note that this last operation 
implies that the qubit is interacting simultaneously with 
the two oscillators. If the interaction is sequential we 
should consider the delay between the two interactions, 
which makes the first oscillator decay for longer (dephas- 
ing can usually be neglected for flying qubits). It is not 
hard to include this effect analytically if needed, however 
our calculations show that it can be taken into account 
with good approximation by doubling the damping rate 
of the oscillators, so that we can present a common treat- 
ment for the two settings. The Wigner function of the 
two oscillators after the qubit is measured in / = e, g is 

W f (a,(3,t,T) = (AA(tT))- 1 (/|^(a,/3,t,T)|/), (36) 

where Af(t, T) = Jd 2 ad 2 (f\W(a, 0,t,T)\f) . 

In the reciprocation procedure the two qubits plus two 
oscillators are described by a 4x4 Wigner matrix, initially 
given by \gg) (gg\ x Wf(a, 0, t, T), and then evolved for a 
further time t (for simplicity the same as before) . Similar 
techniques as above allow to keep the treatment analyti- 
cal. The state of the qubits after reciprocation is obtained 
from the Wigner matrix in both cases of our interest. For 
momentum measurements we use the property 

(p a ,Pb\p\Pa,Pb) = J dx a dx b W(x a + ip a ,x b + ip b ), (37) 

which gives the overlap of the density matrix with the 
momentum eigenstates. In this case we decompose the 
complex variables into their respective real and imagi- 
nary parts, i.e. a = x a + ip a , — x b + ip b . When 
measuring parities we use 

W(a,0) = ^Tr{D a ( a y D b {0)^ pD a (a)D b (0)Tl a fl b }, 

(38) 



which is a well known alternative definition [13j |. For 
example, the two qubit density matrix corresponding to 
the parity eigenvalues II a = Ilf, = +1 is 



p++(t)oc / d 2 ad 2 0W(a,0) + - J d 2 aW(a,0)+ 

/2 
d 2 0W(O,0) + ^W(O,O). (39) 



ENTANGLEMENT ANALYSIS 

We proceed to evaluate the entanglement contained in 
the states prepared at each stage of the protocol. Let 
us start with state pj of the two oscillators after the 
entangling qubit is measured. We investigate its nonlocal 
properties by looking for a Bell inequality violation. We 
consider the measurable quantity 



AP(a,0) = — W f (a,0)=P sl 



3 (a,/3) - Pdiff(a,/3), 
(40) 

where P samc /difi{ot, P) is the probability of measuring the 
two oscillators with same/opposite parities after having 
applied the displacements D a (a) j ', D b (ffi . This quantity 
takes values between -1 and 1 and therefore it can be used 
to form a Bell inequality: 

B(a, P, a', 0') = AP(a, (3) + AP(a', f3)+ 

+ AP(a,0')-AP(a',0')<2, (41) 

where a,0,a',/3' are arbitrary complex numbers. A value 
of B greater than 2 indicates the presence of quantum 
nonlocal correlations. The maximum violation for a bi- 
partite system is B = 2v2, corresponding to strong non- 
locality. In the limit of large t and re, 7 — > 0, the maxi- 
mum value of B is given by 



B u 



1 



(l + 2n(T)y 



:2V2, 



(42) 



corresponding to a set of solutions {a, 0, a', 0'}. We can 
see that if T is lower than a critical temperature T c ~ 
0.408cj, then a violation of the Bell inequality holds. In 
the case of finite times and re, 7 > 0, we get a lower bound 
for the Bell inequality violation by considering 



>B(a,0,a',0'), 



-Kt/2 



(43) 

). By intro- 



where t has to be substituted by -(1- 
ducing decoherence the time window in which it is possi- 
ble to verify a Bell inequality violation becomes limited 
and the maximum value of B drops from the ideal value 
of 2^2- Nevertheless if the decoherence and dephasing 
rates are reasonably small we are still left with a range 
of times and temperatures which give a value of B signif- 
icantly larger than 2, as we can see in Fig. [TJ 

Let us discuss our results for the reciprocation proce- 
dure. In the case of momentum measurements we can 
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a b 




0.4 0.4 

FIG. 1: Lower bound for the Bell inequality violation 
Max(B, 2), relative to state p g as a function of temperature 
and interaction time. The case p e is qualitatively similar. Plot 
(a) corresponds to the ideal case 7, k ~ 0, where the maximum 
violation compatible with the system's temperature is reached 
for long enough times. In plot (b) (k, 7) = (5, 3) x 10~ 3 as in 
0, which limits the time window to circa 2 Rabi periods and 
makes the maximum violation drop to approximately 2.4 



. 




FIG. 2: Entanglement reciprocation via momentum measure- 
ments, for the entangling qubit measured in e. The g case 
is qualitatively similar. In plots (a) and (b) time and tem- 
perature have been fixed (i = 2 and T — uj). Plot (a) shows 
the probability distribution for the outcomes p a ,Pb- Plot (b) 
shows the fidelity of the two qubit density matrix to the Bell 
state |V> + ) = ^(|ff e ) + l e <?))- A fidelity above 0.5 indicates 
an entangled state. Entanglement is conditioned to measuring 
p a ,Pb inside the central peak of the probability distribution, 
which has a volume ~25%. The remaining peaks yield sep- 
arable qubit states. In plot (c) the maximum fidelity (corre- 
sponding to p a ,Pb — 0) is shown as a function of temperature 
and interaction time. We set k — 7 = 0.01 in all plots. 



find, conditionally, the qubits in an entangled state for 
values of T well above T c , as shown in Fig.O Remarkably 
this procedure sustains higher and higher temperatures 
as the parameters k and 7 approach zero. For k, 7 = 
we found no temperature limit at all, and we only had to 
choose long enough interaction times to let the entangle- 
ment appear. Fig[3] shows the entanglement (negativity) 
of the two-qubit density matrix in the case of reciproca- 
tion with parity measurements, averaged over the four 
possible outcomes n a ,Ilb = ±1. In this case we found 
again a limit in the sustainable temperature, which does 
not increase as k, 7 — > 0. The presence of a maximum 
allowed temperature may be a peculiarity of parity mea- 
surements. It is important to point out that the interac- 
tion time does not need to be finely tuned at any step of 
our protocol, provided that we remain in the proper time 
window which allows the entanglement to appear. In all 
the situations presented this time window increases indef- 
initely as we increase our ability to build systems with 



smaller k and 7, so that technological advances in build- 
ing higher quality qubits and oscillators will eventually 
remove any need of fine tuning in the interaction times. 

POSSIBLE IMPLEMENTATIONS AND 
CONCLUSIONS 

Presently cavity-QED setups with flying atoms seem 
to be a promising ground for the implementation of 
our scheme, e.g. in [l[ Rauschenbeutel et al. achieved 
(ft, 7) ~ (0.007,0). Circuit-QED setups with supercon- 
ducting qubits are also very appealing since values of 
(k, 7) — (5,3) x 10~ 3 or even better can be achieved by 
present technology Q. Values of (k, 7) ~ (0.23,0.005), 
which do not allow for a Bell inequality violation but still 
give decent results for reciprocation (e.g a maximum Bell- 
state fidelity above 0.85), seem feasible at present for a 
spin coupled to a nanomechanical resonator, as shown in 
[J] . In the context of quantum nanomechanical systems 
coupled to superconducting qubits k ~ 0.01 is within 
reach, but the rapid qubit dephasing 7 ~ 1 is still an is- 
sue for present technology, although some progress seems 
possible in the near future Note that parity mea- 
surements have already been performed successfully in a 
cavity-QED setting [l4j , while momentum measurements 
should be possible with capacitive transducers for quan- 
tum nanomechanics. 

We are convinced that the robustness of the protocol 
relies on the features of the Hamiltonian ([2|). In fact, as 
we see in ((5J, it allows us to apply to the oscillators a 
coherent superposition of displacements in opposite di- 
rections, with amplitude proportional to t. We expect to 
see entanglement between the two oscillators if t is larger 
than the phase space extension of the initial state, even 
if mixed. However, large t also means more decoherence, 
thus the competition between these two effects will de- 
termine the existence and extension of a time window in 
which entanglement can be established. 

We thank the Engineering and Physical Sciences Re- 
search Council United Kingdom, the Quantum Informa- 
tion Processing Interdisciplinary Research Collaboration, 
the Royal Society and the Wolfson Foundation. 
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FIG. 3: Average negativity of the two-qubit density matrix 
in the case of parity measurements. Plot (a) refers to the 
entangling qubit being measured in g, while plot (b) refers 
to the outcome e. The maxima are close to the value ~ 0.3. 
(k,7) = (0.007,0) as in [l[ 
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APPENDIX: SOLUTION FOR GENERIC INITIAL 
CONDITIONS VIA THE INTEGRAL KERNELS 

The integral Kernels Vj(a,ao,i); j = 1, ..,4; are solu- 
tions of the equations (|16H19p with initial condition 



Vj(a, ao, 0) = S(a — a Q ). 



(A.l) 



From the integral kernel the solution corresponding to 
an arbitrary initial condition Vj(a, 0) can be expressed in 
terms of a convolution integral: 



Vj(a,t)= / Vj(a,ao,t)vj(ao,0)da 



(A.2) 



It is possible to derive the four integral Kernels as fol- 
lows. We know from the literature (see e.g. 13]) that 
the function 



K (a,a ,t) 



1 



ttA(1 



') 



A(l- 



(A.3) 



is the integral kernel for L(T), i.e. 



d t K (a,a ,t) = L(T)K (a,a ,t), (A.4) 
K o (a,a o ,0) = S(a- a ) (A. 5) 



To find the four integral kernels of equations (|16H19j) , we 
consider the four ansatz 

Vi(a,a ,t) = K (a + \(t),a ,t) (A.6) 
V 2 (a,a ,t)=Ko{a-\{t),a ,t) (A.7) 
V 3 (a, a ,t) = K (a, ao ,t)e^ a+ao+a " +a ° )+v ^ (A.8) 
Vi(a, a ,t) = K (a, a(ht ) e ~^ a+ao+a ' + <^ w 

(A.9) 



Using the techniques previously described, we see that 
the first two ansatz lead again to equation (|25p . which 
we already solved. For the last two ansatz to be ac- 
tual solutions, the following differential equations must 
be satisfied: 



2i 



fl, = K 



K 1 H 

2 T~- 

e -Kt/2 



V 



1 -e 

2 



— Kt 



nAfi — 7. 



(A.10) 

(All) 
(A.12) 



These are obtained by substituting (|A.8|) and (|A.9|) into 
(HHJ) and (Unj) respectively, then collecting separately the 
terms proportional to Vj,aVj and a^Vj, with j = 3,4. 
The system of equations above is overdetermined in gen- 
eral, however it is possible to find a particular solution 
which also accounts for the initial conditions fj,(0) = 
i/(0) = 0: 



/i(i) = 



4i 1 



,-Kt/2 



k 1 + e 



-Kt/2 ' 



/(<) = kA 

Jo 



/x(r) dr — 7*. 



(A.13) 



(A.14) 



We have at this point all the tools needed to describe the 
evolution of any given initial state in terms of integrals. 
The procedure is easily generalizable (at least in princi- 
ple) to an ensemble Q = {q 1 , q^}, containing k qubits, 
and an ensemble A — {a±, a{\ containing I oscillators, 
such that each qubit q,j is coupled to a subset of A with 
the Hamiltonian @, while objects of the same kind do 
not interact. 
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